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Abstract. One way to account for the acceleration of the universe is to modify 
general relativity, rather than introducing dark energy. Typically, such modifications 
introduce new degrees of freedom. It is interesting to consider models with no new 
degrees of freedom, but with a modified dependence on the conventional energy- 
momentum tensor; the Palatini formulation of f(R) theories is one example. Such 
theories offer an interesting testing ground for investigations of cosmological modified 
gravity. In this paper we study the evolution of structure in these "modified-source 
gravity" theories. In the linear regime, density perturbations exhibit scale dependent 
runaway growth at late times and, in particular, a mode of a given wavenumber goes 
nonlinear at a higher redshift than in the standard ACDM model. We discuss the 
implications of this behavior and why there are reasons to expect that the growth 
will be cut off in the nonlinear regime. Assuming that this holds in a full nonlinear 
analysis, we briefly describe how upcoming measurements may probe the differences 
between the modified theory and the standard ACDM model. 
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1. Introduction 

The concordance cosmological model describes a universe consisting of approximately 
5% ordinary matter, 25% dark matter, and 70% dark energy. While it fits a wide variety 
of data, this model relies heavily on the existence of a "dark sector" comprising 95% 
of the universe. Given the mysterious nature of dark matter and dark energy, and the 
fact that their existence is inferred exclusively through their gravitational effects, it is 
natural to wonder whether the apparent need for these components could be a sign that 
gravity is deviating from conventional general relativity (GR) on large scales. 

Dynamical measurements that are taken to imply the existence of dark matter 
generally refer to length scales of kiloparsecs or greater, while evidence for dark energy 
comes from the acceleration of the universe, a phenomenon characteristic of the present 
Hubble radius of order one gigaparsec. The most precise experimental tests of GR, 
meanwhile, probe much smaller length scales [Tj; Solar System measurements cover less 
than a milliparsec, while the binary pulsar PSR 1913+16 has an orbital semi-major 
axis less than a microparsec. We can therefore imagine that there is a large dynamical 
range over which gravity obeys Einstein's equation to high precision, while behaving 
differently at sufficiently large scales. 

Within the modified-gravity category there have been several different proposals: 
induced gravity in a class of extra dimensional models (DGP braneworlds) [21 EH Ej; 
phenomenological modifications of the Friedmann equation of cosmology [21 E]; and 
direct, manifestly covariant modifications of the four-dimensional Einstein-Hilbert 
action EDI [HI 112 113 El 113 • (For some other ideas see HBJ HH HH1 EEH] .) In the 
last of these approaches, the simplest models involve an action obtained by integrating 
a function of the curvature scalar R, 



and varying with respect to the metric g^ v . These models were soon discovered to 
be dual to scalar-tensor theories, which lead to unacceptable deviations from GR in 
the solar system (201 ■ However, certain models based on inverse powers of more general 
curvature invariants (particularly those involving the square of the Riemann tensor) 
have recently been shown to agree with solar system tests [21] and to provide a good fit 
to the supernovae data [22] > although tight constraints arise from the requirement that 
they be ghost-free [23 | l2l l 125] . 

On theoretical grounds, however, there are obstacles to constructing infrared 
modifications of GR that will escape detection on smaller scales. In a weak-field 
expansion around flat spacetime, GR describes the propagation of massless spin-2 
gravitons coupled to the energy- momentum tensor T^, including that of the gravitons 
themselves. But such a theory is essentially unique; it has long been known that we can 
start with a field theory describing a transverse-traceless symmetric two-index tensor 
propagating in Minkowski space, couple it to the energy-momentum tensor of itself and 
other fields, and demonstrate iteratively that the background metric disappears, leaving 




(1) 
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instead a curved metric obeying Einstein's equation. It is therefore generally believed 
that infrared modifications of GR will necessarily involve the introduction of new degrees 
of freedom. 

For the purpose of explaining the acceleration of the universe, however, there is 
a loophole in this argument. The properties of gravitons define both the propagation 
of linearized gravitational waves and the Coulomb form of the Newtonian gravitational 
potential in the weak- field limit. But there is more to gravity than gravitons, even in 
the infrared. In particular, the evolution of a Robertson- Walker cosmology is described 
by the Friedmann equation, 

87rG k 
3 a 1 

where p is the energy density, k is the spatial curvature, and H is the Hubble parameter, 
related to the cosmological scale factor a(t) by H = a/ a. The Friedmann equation has 
nothing to do with gravitons; it is a constraint, relating the instantaneous expansion 
rate to the energy density and curvature. In fact, the Friedmann equation is a particular 
example of the Hamiltonian constraint of general relativity, relating the embedding of a 
partial Cauchy surface in spacetime to its intrinsic geometry and the energy-momentum 
tensor. 

Because the Friedmann equation arises as a consequence of the constraint structure 
of GR rather than from the behavior of gravitons, we are free to imagine modifying it 
without introducing new degrees of freedom. In particular, Einstein's equation relates 
the Einstein tensor to the energy-momentum tensor for matter, 

Gfiv = StcGT^ u . (3) 

In a theory describing a set of matter fields x% with an action S^lg^, x%], the energy- 
momentum tensor is defined by 

T^ S -4^. (4, 



5g^ 

We may ask, then, to what other sources could gravitons couple? In other words, we 
might define an alternative theory obeying 

G^ = 8nGT^(xi) , (5) 

where is some symmetric (0, 2) tensor that is conserved by virtue of the matter 
equations of motion, but differs from the conventional form (J3J). Since we have altered the 
right-hand side of Einstein's equation without introducing any new degrees of freedom, 
it may be possible to account for the acceleration of the universe without dark energy 
while remaining consistent with conventional tests of GR. 

In fact this idea has been realized in the form of the Palatini formulation of f(R) 
gravity. As Flanagan has shown JU] 5 the equations of motion obtained by separately 
varying the metric and connection in an action of the form actually leads to a theory 
with fewer degrees of freedom; there is no propagating scalar, only the massless spin-two 
graviton of ordinary GR. In this paper we further consider models of this type, leaving 
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behind the inspiration of the Palatini formulation of f(R) models and working directly 
with theories that have no propagating scalar, which we dub "modified-source gravity" 
(MSG). The theories we consider are not precisely identical to the Palatini models, 
since we do not separately vary the metric and connection, so that the equations for 
the matter sector will be slightly different; nevertheless, there is an essential similarity 
between the two ideas. 

Flanagan has argued that models of this type are experimentally excluded, as 
the effective matter action includes higher-order couplings not seen in particle-physics 
experiments. However, as we argue below, there are a number of reasons to believe 
that these constraints can be avoided; furthermore, the MSG field equations are an 
interesting toy model of cosmological modified gravity against which observations may 
be compared. 

The possibility of distinguishing between modified gravity and dark energy by 
comparing the expansion history of the universe to the growth rate of cosmological 
perturbations has recently been emphasized from a variety of approaches [2E1 EH EHl 
E0lEIlElE3ElESlEniEIlEHlEi- Given the paucity of viable alternative theories of 
cosmological gravity, it is important to develop a model-independent perspective on the 
ways in which this distinction can manifest itself; we know of no better way to develop 
such intuition than to examine as many models as possible. 

In this paper we therefore explore the cosmology of modified-source gravity by 
studying cosmological perturbation theory and extracting predictions for the growth 
of large-scale structure. We find that there is a substantial boost in the growth of the 
gravitational potential in comparison with ordinary ACDM. We discuss the implications 
of this behavior and why there are reasons to expect that the growth will cease in the 
nonlinear regime. Assuming that this holds in a full nonlinear analysis, we briefly 
describe how upcoming measurements may probe the differences between the modified 
theory and the standard ACDM model. 

Modified-source gravity is examined very much in the spirit of a toy model, as a 
laboratory to help understand the possible differences between theories of dark energy 
and theories that alter general relativity. It does not attempt to solve the cosmological 
constant problem (we set the vacuum energy to zero by hand), nor does it require 
less fine-tuning than any other typical model of dark energy (we choose the potential 
delicately so that acceleration happens at very late times), or is there any obvious way 
in which it would naturally appear as the low-energy limit of a more complete theory in 
the ultraviolet. Nevertheless, it is sufficiently difficult to find cosmological alternatives 
to general relativity that a simple explicit model can be useful in helping to focus efforts 
to distinguish between modified gravity and sources of dynamical dark energy. 
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2. Modified-Source Gravity 

2.1. Scalar-tensor and f(R) theories 

The Einstein-Hilbert action for general relativity is 



(6) 



where R is the curvature scalar and Mp = 1/ v8ttG is the reduced Planck mass. Matter 
fields couple minimally to the metric g^. The propagating degrees of freedom of this 
theory are the two polarization states of a massless spin-two graviton. 

One of the simplest ways to modify GR without introducing new fields is to consider 
actions that depend non-linearly on R, 

S f = J d'x^f(R). (7) 

A particular example in which deviations from GR become important at small curvatures 
is f(R) = (Mp/2)(R — fi 4 / R) [3IHj. Interestingly, however, the linear Einstein-Hilbert 
term (0) is the only function of R that propagates a spin-2 graviton by itself; any other 
function f(R) gives rise to a scalar-tensor theory, with both a spin-2 graviton and a 
spin-0 scalar [¥UJ El EE]- This can be seen explicitly by considering an action with 
gravity coupled to a scalar A, 

S = J d 4 xv^ [(it! - A)f (A) + /(A)] , (8) 

where /'(A) = df/dX. The field A functions as a Lagrange multiplier, whose equation of 
motion sets A = R. Plugging this back into (jSJ) yields (J7|), so long as /" ^ (that is, if 
/ is anything other than linear). Defining a new scalar field if) via 

4> = \WW], (9) 

we can perform a conformal transformation of the form 

^„ = e 2V V. (10) 

The conformally-transformed action becomes that of ordinary general relativity coupled 
to a propagating scalar field, 



S = d 4 x\/— # 



M 2 ~ ~ 

-f R - 3sT(W0 a - uty) 



where (V^) 2 = g^V ^ipV u il) . This is the Einstein-frame version of the model, in which 
the Ricci curvature appears by itself, but matter fields couple to the scalar ip through 
the combination g^ v = e~ 2 ^g^ v \ in the matter frame given by (J7J), there is no direct 
coupling between if) and the matter fields, and free particles move along geodesies of 
g^y. The potential is given by 

gW = A/ ' ( 2 A ^ A) /(A) M|. (12) 
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In this expression, A is taken to be a function of if> via inverting (JUJ). (Note that our 
notation differs from that in [7j.) 

The problem with such a model is that, if we choose the function f(R) so that 
gravity is modified at small curvatures and the universe accelerates at late times, the 
scalar field will be very light, and the theory comes into conflict with Solar-System 
tests of gravity. In terms of the Brans-Dicke parameter u, the model of [7j is nearly 
equivalent (apart from the small and presumably negligible potential) to a theory with 
uj = 0; meanwhile, the best current limits come from measurements of the Shapiro time 
delay from the Cassini mission, and give u > 40,000 |43j . However, recently some 
authors have claimed that the GR limit of the dynamical equivalence between f(R) and 
scalar-tensor theories is not well behaved, since f"(R) — > 0, and some f(R) theories 
behave in the Solar System in a manner perfectly consistent with current experimental 
limits IH US EE!- 

Different strategies for avoiding this constraint have been explored. One approach 
is to consider models with inverse powers of other curvature invariants such as R^ V R^ V 
and R^ pa R^ pa [13 I2JJ 122 123 121 123 • As mentioned in the introduction, the degrees of 
freedom of such models differ from those of simple f(R) theories, and may be consistent 
with solar system tests of gravity. Alternatively, one may consider the Palatini versions 
of f(R) theories, in which the metric and connection are treated as independent variables 
|3 IHH EEH U21 HI| ■ Unlike in the case of the Einstein-Hilbert action, f(R) theories do 
not have identical equations of motion in the Palatini formulation as in the conventional 
scenario based on the metric alone. Indeed, Flanagan has shown that the scalar degree 
of freedom disappears entirely [TO] - In the next section we explore a possibility of this 
form, not through the Palatini variation of an f(R) action, but simply by eliminating 
the scalar kinetic term by hand (although of course the formulations are related). 

2.2. Eliminating the scalar 

Our goal is to take the scalar-tensor action (jllj) . equivalent to (|7|) . and eliminate the 
propagating scalar degree of freedom. We choose the most brutally direct approach: 
simply erasing the kinetic term Sg^V ^ipV from We are left with a new theory 
in which tp is a Lagrange multiplier, without any dynamics of its own. Alternatively, 
we could imagine multiplying the kinetic term by a constant, and then taking the limit 
as the constant went to zero; in that limit, the scalar is still propagating, but decouples 
from any other fields. Generally, radiative corrections would tend to drive this constant 
to order unity. We will nevertheless set it to zero for purposes of this paper, accepting 
this as one of the fine-tunings that inevitably accompanies models of dynamical dark 
energy .| 

The full Einstein- frame action for the theory, including additional matter fields Xii 

% Of course, the kinetic term may always be canonically normalized by a field redefinition. Setting the 
coefficient of the kinetic term to zero is therefore equivalent to taking the potential and the couplings 
all simultaneously to infinity. 
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is 



S 



-g 



(13) 



In the absence of the matter action, this model would be completely trivial. The scalar is 
non-propagating, and its equation of motion would fix ip at any allowed value ipo at which 
dU /dip = 0. The model is then simply general relativity with a vacuum energy given 
by U(ip ). The coupling to matter, however, leads to important consequences. We can 
undo the conformal transformation to return to the matter-frame metric = e -2 ^^,, 
yielding 



S 



d 4 x, 



M 



^R+3e^(Vip) 2 -U(iP) 



S(m)[9fiu, Xi] 



(14) 



(15) 



where the new potential is 

U{$) = e^U(ip) . 

The actions (fT4*|) or f)13|) define modified-source gravity. We have left behind our 
original inspiration from f(R) gravity, so the potential U(ip) is now simply a free function 
that defines the theory. In the form (|14p. the model resembles a conventional scalar- 
tensor theory, with a kinetic term for ip and a direct coupling to the curvature scalar. 
However, there are implicitly derivatives of ip in the |e 2 ^i? term, which would be revealed 
after integration by parts, and would cancel the explicit kinetic term. The scalar field ip 
is actually completely non-dynamical, as is evident from the Einstein- frame expression 
(|13|). in which no derivatives of ip appear. In principle, this field could be integrated 
out exactly, and we will examine this approach in the next section; for the moment, it 
is more convenient to leave ip explicitly in the action and field equations. 

Let's examine the theory in the matter frame ()14j) . The gravity equation of motion, 
obtained by varying with respect to g^, can be written 



Here, T, 



WO 



fxu is the effective energy-momentum tensor for ip, 

- [U(iP) + (V^) 2 + 2Uip] 9ia , - 2V^V u iP + 2V fl V u iP 
The equation of motion for ip is 
' _, --< 

dip 



where □ = g^V '^V '„ 

Dip+ (Vlp) 2 + 



6M 2 



*™ - 1 -r 

6 







(16) 



(17) 



(18) 



Again, these equations bear a close resemblance to those of ordinary scalar-tensor 
gravity. However, we can take the trace of ()16|) to obtain 

R = ~mJ [ ' t + Aum + 6(w)2 + ^ ' (19) 

where T = g^ u T^ is the trace of the matter energy-momentum tensor alone. 
Subtracting this from the scalar equation (fT%|) leaves 
dU 

cr0 



AU{ip) 



-T . 



(20) 
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Thus, the trace of the gravity equation can be used to eliminate all spacetime derivatives 
from the scalar equation, leaving us with a purely algebraic equation for ip in terms of 
the matter fields. 

If we start with U (ip), we could use this equation to get T(ip), and invert that to get 
if)(T). If we start with ip(T) [or equivalently T(ip)], we can then express the potential 
in terms of ip as 

U{$) = -e 4 ^ / e-^'T{i>') dip . (21) 

Jtpo 

Lastly, if we know U(T) instead of ip{T), we could get ip via 

V>CO = - Jf(j& + 4U )* r - ( 22 ) 



3. Comparison with Experiment 

3.1. Solar-System tests 

Consider the vacuum equations, where T^ v = 0. Then from pf)|) the scalar is pinned at 
some value ipo, with 

4<7(^ ) - U'M = . (23) 

Looking back at the vacuum gravitational field equation is 

MpG^ = -e-^U(ip ) gilu . (24) 

Thus, e~ 2 ^°U (ipo) plays the role of an ordinary cosmological constant. Otherwise, the 
vacuum equation is precisely the vacuum Einstein equation. The Schwarzschild metric 
is therefore an exact solution of this theory (if the vacuum energy vanishes; otherwise 
it would be Schwarzschild-de Sitter). Gravitational waves in free space are completely 
identical to those in general relativity; the propagating degrees of freedom are simply 
those of a massless spin-2 graviton. 

In f(R) theories, Schwarzschild (for which R^ v = 0) is generally also an exact 
solution to the vacuum field equations, even though those equations are not identical 
to Einstein's. However, it is not a unique solution, even with spherical symmetry; 
Birkhoff 's theorem, which relies on Einstein's equation, fails to apply. This failure can 
ultimately be traced to the existence of new degrees of freedom; gravitating bodies 
are not only characterized by their mass, but also by a charge that couples to a new 
scalar field. Modified- source gravity, in contrast, has no new degrees of freedom, and 
Einstein's equation holds in vacuum, so Birkhoff's theorem applies. The Solar-System 
tests of gravity that rule out simple f(R) theories are completely compatible with MSG. 

We are free to choose the "potential" U(ip) to satisfy phenomenological 
requirements. If our interest is in causing the universe to accelerate at late times without 
affecting local tests of gravity, these include the following: 
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• For energy densities p larger than some critical value p* of order the present mass 
density of the universe, we want to recover conventional cosmological expansion 
(a oc t 2 l n when the total energy density is evolving as p oc a~ n ). This is achieved 
if ip — > ip* = constant at p ^ p*. For constant ip, the additional energy- momentum 
contributions from (fTTj) are simply those of a constant effective vacuum energy 
U{ip) < p*. In particular, the effective gravitational constant will be constant, and 
we recover conventional Friedmann-Robertson- Walker cosmology 

• Since p* is smaller than the density of any compact astrophysical object, the 
previous requirement guarantees that we will automatically get xjj = ip* inside 
planets and stars. The strength of gravity will not be significantly density- 
dependent, so the predictions of Solar-System tests or the binary pulsar are 
unaffected without further restrictions. 

• For p ^ p* the universe should accelerate. In section 14.21 we consider one concrete 
choice of U(ijj) that accomplishes this goal. We might also ask that the value of U 
vanish when p = 0, so that Minkowski space is a solution to the theory; otherwise, 
late-time acceleration could simply be blamed on the new potential energy, which 
wouldn't really be different from ordinary dark energy. 

3.2. Matter interactions 

Despite the characterization of MSG as a theory of "gravity," we can certainly imagine 
integrating out the scalar field in the Einstein-frame action (jlHj) to obtain a modified 
matter Lagrangian coupled to general relativity. The result will be an action written 
purely in terms of the ordinary matter fields and gravity; if the matter action looks 
conventional in the matter frame, it will generally be an ungainly mess in the Einstein 
frame with ip integrated out. As noted by Flanagan ^Uj, this implies the existence of 
higher-order (nonrenormalizable) terms in the matter action, which will be proportional 
to inverse powers of a very small mass scale p characterizing the regime in which 
gravity is to be modified. Since p should be very low in models relevant to late-time 
acceleration, such new interactions would presumably be easily noticeable in experiments 
(and, needless to say, have not been seen). 

In addition to the introduction of new nonrenormalizable interactions, another 
phenomenological problem identified by Flanagan jTUj relates to the nonlinearity of the 
source for gravity in MSG. The effective energy density is a nonlinear function of the 
conventional matter-frame energy-momentum tensor. Thus, it is illegitimate to average 
the energy density of a particulate medium over some coarse-graining to obtain a fluid 
description. In the limit where a body is made of pointlike particles, the density is 
infinite at the locations of the particles and zero in between, and treating it as a smooth 
density profile would be a mistake. 

Despite these reasonable concerns, for the rest of this paper we will imagine that 
MSG is experimentally viable and an averaged fluid description for the matter fields is 
sufficient. Our primary justification for sidestepping this important issue is the dramatic 
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separation of scales between the regimes we are considering (evolution of cosmological 
perturbations) and those in which the above considerations become relevant (atomic and 
particle-physics scales). We therefore find it interesting to consider the MSG equations 
as a phenomenological description of gravity in the infrared, even if we do not have a 
viable model on small scales. In order to best understand how we can observationally 
distinguish dark energy from modified gravity, it is important to characterize different 
ways in which modified gravity could manifest itself cosmologically. Since the MSG 
equations provide a consistent dynamical description of gravity on very large scales, 
studying their observational consequences is a useful practical exercise. 

Furthermore, it is certainly possible to imagine ways in which the problems 
of higher-order interactions and the non-linear energy- momentum tensor could be 
overcome. Perhaps the most straightforward would be to imagine that the kinetic term 
for our auxiliary scalar if> were not exactly zero, but rather some very small number; 
the fluctuations of this field could then serve to average over the density of particles on 
sufficiently small scales. Alternatively, if the dark matter is some smoothly-distributed 
bosonic condensate (such as axions), the fluid description used in this paper could be 
completely accurate. In fact, it is important to note that experimental results relevant to 
the question of higher-order interactions do not actually take place in a true vacuum, but 
rather in the presence of whatever dark matter background exists in the Solar System. 
Since the local density of dark matter is likely to be substantially greater than that 
of the average density of the universe, terrestrial experiments take place in a regime 
where gravity is accurately described by GR and novel MSG effects are suppressed. 
Finally, Flanagan also notes jTTj that the new interactions are strongly-coupled in the 
low-energy regime in which they are purportedly observable. All of these possibilities 
are interesting and subtle, and deserve greater attention than we will give them in this 
paper, where our interest is exclusively cosmological; further development of the theory 
to put it on more secure microphysical foundations would be very worthwhile. 

4. Robertson- Walker Cosmology 

4-1- The Modified Friedmann equation 

In this subsection we examine the evolution of a homogeneous and isotropic universe, 
deriving the modified Friedmann equation for MSG. In the next subsection we focus on 
an explicit choice for the potential U (if)) that leads to late-time acceleration. 
Consider a flat Robertson- Walker metric, 



where k is zero, positive or negative depending on the curvature of the spatial 
hypersurface. It is not normalized, as we have chosen instead to set a = 1 at the 
present time. The matter fields are taken to be a perfect fluid, with energy-momentum 




(25) 
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tensor T^ v = di&g(—p,p,p,p). The 00 component of the gravitational field equation 
(fTSJl becomes 

-2'tp 



3H " + ^ = Z mJ [p + um ~ 3 ^ " 6 ^ ' (26) 

where H = a/a. We can convert the time derivatives of if) into derivatives with respect 
to a by using 

^ = M = h(^-). (27) 



(28) 



da \ d In a 

Collecting terms proportional to if 2 and dividing, we obtain 



H 2 



dip 



-2V- K 



din a, 

which serves as the cosmological evolution equation for MSG. 

There are two obvious modifications from the conventional Friedmann equation: a 
potential-energy contribution U(ip), and a variable-strength effective Newton's constant, 

^SK 1 ^)"- (29) 

When re = and p ^> U, this is the quantity that relates the energy density to the 
expansion rate. In fact, both the potential and the effective Newton's constant are 
simply functions of the density p and the pressure p through the if) equation (|20jl. The 
new equation is therefore reminiscent of the Cardassian model J3J, in which the right- 
hand side of the Friedmann equation is a non-linear function of p. 

4-2. A model 

There is a great amount of freedom in the choice of MSG dynamics, represented by the 
arbitrary function U(ip). In this subsection we discuss one simple example that satisfies 
the criteria laid out in section I3.ll including late-time acceleration. 

For simplicity we imagine that the matter is represented by a pressureless fluid, so 
that if) can be thought of as a function of the matter energy density 

p = -T = p a- 3 , (30) 

where po is the average cosmological energy density today. Given the relations ([20)1 . 
(J2T|) . and (J22J) . we can choose to specify any of the three quantities ip, U, and p as a 
function of any one of the others, and their relationship to the third will be automatically 
determined. 

We would like to choose behavior for which if) is approximately constant when the 
energy density is substantially higher than the present average cosmological density po, 
and decreases for p p Q . A convenient form to choose is 

-iip _ ^ ( Po 



e-^ = a(^H 1. (31) 
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with a a dimensionless parameter defining the model. The limiting behavior is given by 

1 



ip(p 



oo 







V>(p-0) = -lii(p) 



-oo, 



(32) 



satisfying the aforementioned criteria. 

The convenience of this model arises from our ability to analytically determine the 
potential as a function of the density, 



TTl \ a P°P 1 ( P° 

U ip) = ——, In a — 

4(ap + p) V P 

or equivalently clS 8b function of ip, 



As p — > we get 



^- i ln(l-e 4 ^) 



U(p^0) = P -\n(P- 
4 VPo 



0. 



(33) 



(34) 



(35) 



The fact that the potential vanishes at zero density guarantees that Minkowski space is 
a solution to the model. At large density we have 



U(p 



oo 



(36) 



4 VPo. 

which rises more slowly than p, so that the potential does not dominate at early times. 
The behavior of U and ip as functions of p is shown in Figure ^ 



a. 




10 10 




Figure 1. The potential U (in units of po) and scalar field ip (in units of Mp) as 
functions of the cosmological energy density p. 



The other factor appearing in the Friedmann equation (|28j) involves dip /din a, which 
can be written in this model as 

dip dip 3ap 
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Putting everything together, the full Friedmann equation (J28|) for the specific model 
defined by (}3~T]) becomes 



H 2 



4ap + 4p 
ap + 4p 



3iW| 



ap + p - (apo/4) ln(ap /p) 



^2 



(3* 



It seems unlikely that such an expression would have been guessed at as a 
phenomenological starting point for exploring cosmological dynamics in theories of 
modified gravity. 

An observer, understandably assuming the validity of conventional general 
relativity, would interpret measurements of the expansion history a(t) in MSG in terms 
of the derived density and dynamics of an effective dark-energy component p^ E . In terms 
of the (in principle) observable quantities H and p, the effective dark-energy density is 

f& = 3M^ 2 - p . (39) 

where the matter density of course evolves as p oc a~ 3 . The effective equation-of- state 
parameter is 

l ddnpf| 
3 din a 



■1 



(40) 



The behavior of these quantities as a function of redshift in the model defined by 1)31)) 
is shown in Figure El 
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Figure 2. The plot on the left shows the effective dark-energy density, in units 
of the matter density today, that observers would reconstruct. (That is, the dark- 
energy density that would lead to equivalent behavior for the scale factor if general 
relativity were correct.) The right plot shows the corresponding effective cquation-of- 
state parameter as a function of redshift. 



We fix the parameters of our model by fitting the luminosity distance-redshift 
relation to the Supernovae Legacy Survey (SNLS) data set of 115 type-la supernovae 
|49j : in addition we fit for the distance to the surface of last scattering: we use the 
third-year WMAP results jSO] to fix its redshift, z\ ss = 1088^) the acoustic peak scale, 
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I a = 302^4, and its calibration through the matter density, Vt m h 2 = 0.1265 ± 0.008. 
In doing the fit, we include both matter and radiation, with a total energy density 
P = Pm (l + — ), where R is the radiation to matter energy density ratio today. In 
fitting to the data, we allowed for the existence of spatial curvature, characterized by 
the parameter 



C 



rm 



(41) 



where is the matter energy density today. The fit prefers a slightly open universe, 
but achieves approximately the same x 2 as ACDM. The set of best fit parameters is 



a 



0.98 



C = 0.12 



h = 0.72 



(42) 



Since we are fitting cosmological parameters in our model to the data, we need a 



consistent definition of the fractional densities Q m and fi„ 
defining the critical density today in the usual way, via p° a 
write the Friedmann equation (|38j) as 



To do this, we begin by 
- 3Mp\Hq, enabling us to 



H 2 



16(cra 3 + l) 1 



(era 3 + 4) 2 

Now consider the asymptotic form of this equation as a 
logarithm, we obtain 



4 (1 + f ) (era 3 + 1) - era 3 ln(cra 3 
4a 3 Vera 3 + 1 



C 

+ ^7 



^.(43) 

rev 



H 2 = H 2 



, R \ 1 C 

1 + - ^ + ^7 



a 1 a° 



0. Ignoring the slowly- varying 



(44) 



which then allows us to consistently define the fractional densities today in precisely the 
same way as in ACDM: 

c P v ; 



a 



pI 

r cr 



0.26 



a. 



p° 
r cr 



+0.02, 



(45) 



where the numbers are the best-fit values to the supernova and CMB data. 



5. Linear Perturbations 



To study perturbations, we introduce small time- and space-dependent deviations from 
the background cosmological solutions for the metric and matter sources. Applying 
the equations of motion and linearizing, we then obtain a set of coupled first-order 
differential equations describing the coupled evolution of matter, radiation and metric 
perturbations as the universe expands. 

The evolution equations for matter and radiation in the presence of a perturbed 
metric are given by the Boltzmann equations (for a particularly clear description of this 
see [48]). This formalism is independent of the equations of motion for the metric itself; 
moreover, as we work in the matter frame, the Boltzmann equations are formally the 
same as those derived in standard general relativity. However, these equations do contain 
new dynamics even if their structure is unchanged, since the background quantities on 
which they depend are solutions of modified background equations. 
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From (|42|) . the best-fit background cosmology requires a small negative curvature. 
However, this will negligibly affect perturbations. Therefore, although we may include 
curvature in the background evolution, we neglect it in the treatment of perturbations. 
We consider scalar perturbations to a flat FRW metric in conformal Newtonian gauge; 
the perturbed line element can be written as 

ds 2 = -(1 + 2*(f , t))dt 2 + a 2 (l + 2$(f, t))dx 2 , (46) 

where <£> and \I/ are spacetime-dependent gravitational potentials. We assume that the 
universe contains both radiation and dark matter, which allows us to parameterize the 
perturbations to the energy momentum tensor by 

T° = -p5 (47) 
T°i = (1 + w)pd iq (48) 
T\ = wpi&j + ttV) . (49) 

Here, p is the background energy density, 5 = Sp/p, w is the equation-of-state parameter 
for the background fluid, q is the momentum density, and n l j is the anisotropic stress. 

We can then linearize the modified Einstein equations (|T6*|) and obtain the Poisson 
and anisotropy equations for the perturbed cosmology in Fourier space, 

+ 5ip) + 3 [H + ^) ($ + 5ip) - 3 (H + iy * 

= 2M^( p6 + {U ' ~ 2U ~ 2P)5 ^ (50) 
a 2 e -2 ^ 

$ + VI/ = ~P2M!; w/wr " 2 ^' (51) 

where ip and U are the background quantities and Si/j is the perturbation to if}. The 
scalar anisotropy 7r is related to the anisotropic stress by 7i l j = (d l dj — |5]5 2 ) it, and is 
non-zero only for radiation. We therefore neglect n and from now on. 

Note that, in comparison with the unmodified (ACDM) version, the anisotropy 
equation (|51j) contains just one extra term, namely that involving 5if>. However, this term 
plays a crucial role in defining the difference between dark energy and modified gravity. 
Bertschinger [22] has recently emphasized that knowledge of both the modification to the 
Friedmann equation and the modification to the evolution of the anisotropy, is key 

for the understanding of the effects of the new gravity theory on cosmological structure 
formation. It is the latter which differentiates modified-gravity theories from models 
with dark energy, and therefore it is the existence of this new term which ultimately 
allows us to test the modified-gravity idea. 

In contrast with the anisotropy equation, the evolution equation ()50|) has several 
new entries arising from the modification to the Einstein equation f)lfij) . Besides a 
rescaling of Newton's constant, there are terms that are functions solely of the scalar 
field, its perturbation and the potential. However, by using equation (j2fl|l . we can 
express if> and Sip in terms of the matter density and its perturbation. In particular, we 



Modified- Source Gravity and Cosmological Structure Formation 



16 



have 



Sip 



1 



(52) 



3 d In a 

where 5 m is just the matter density contrast, since radiation makes no contribution to 
T = —p + 3p. 

The time-space (0,j) component of the linearized Einstein equation is 

e" 2 ^ 

$ - Hm = 77772"P(1 + w)q + ipty - Sip + H5tp + ipSip . (53) 

Combining this with ()50|) and making use of (|52|). we obtain an algebraic equation for 
the potential At late times, we are interested in modes well within the Hubble radius, 
and we can also neglect contributions from radiation. We then have 

/ dip \ _ k 2 dip 



a 2 



2M 2 ! 



1 + 



Pu 



d In a J r 111 3a 2 d In a 
This constraint relates the potential \I> directly to the matter variables. 



(54) 



6. The Growth of Linear Structure 



In the previous section we have derived the basic linearized equations for the evolution 
of perturbations in the context of modified- source gravity. In this section we study the 
growth of linear structure. 

At late times, when all the modes of interest have entered the horizon, and radiation 
and momentum flow are negligible, we can combine (|54j) with the Boltzmann equations 
for dark matter to obtain a second-order differential equation governing dark matter 
perturbations 

- 2r / dip \ _ k 2 dip 



S m + 2H8 m - 



2M 2 l 



1 + 



d In a 



3a 2 d In a 



s„ 







(55) 



The first term in the bracket multiplying 5 m is the same as the term in the Poisson 
equation in GR, modified by the fact that the value of Newton's constant is evolving 
as the average density in the universe decreases beneath the critical value at which the 
modifications to GR become important. 

The second term in the coefficient of 5 m in (|55jl is very different in nature: it 
introduces a scale dependence in the growth of structure. For negative dip /din a, small 
scales will begin to grow more quickly than large scales once the universe approaches 
the accelerating phase. In figure El we give an example of this enhanced growth for 
the choice of potential corresponding to the model defined by (j31j) with the best fit 
parameters (|32~j) . We have plotted the exponent n of 



5 m = 5 exp 



din a n(a) 



(56) 



with 5q cl constant. The rate of growth is strongly scale-dependent once the universe 
enters the epoch of acceleration, with the exponent departing from n — 1 (holding for all 
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Figure 3. The exponent of the growth of matter density perturbations, as defined 
in ((56(1 for a selection of modes, solved for the model defined by 1)31(1 with the best fit 
parameters 1(42(1 f near-horizon effects are ignored). The rate of growth is very strongly 
scale-dependent once the cosmology enters the epoch of acceleration, with the exponent 
departing sharply from n — 1 and tending to n oc k at late times. 



modes during matter domination) and tending to n oc k at late times. The fc-dependence 
of the growth rate will result in smaller scales reaching non-linearity extremely quickly 
after the onset of acceleration. Thus, we generically expect that the matter power 
spectrum will be non-linear at larger scales than in ACDM. 

In fact, this enhanced growth of scales which are small but still in the linear regime 
seems to be a generic feature of modified-source gravity. The effect can be traced to the 
existence of the new fc-dependent term in equation ()56|) for 5 m : small scales will witness 
enhanced growth so long as the potential U{ip) is chosen so that d^/dlna < 0. But, as 
inspection of shows, it is precisely this behavior that makes the universe accelerate 
at late times, by increasing the effective value of Newton's constant in the modified 
Friedmann equation. It therefore seems difficult to avoid this phenomenon simply by a 
clever choice of the potential U (ip) . 

A possible loophole in this argument would be to consider models in which ip were 
nearly constant in the present era, but the contribution of the potential U(ip) itself 
to the right-hand side of were to induce cosmological acceleration. One might 
reasonably complain that this case is simply a dark-energy model in disguise, rather 
than a modification of gravity; however, we should keep in mind that U(tp) is not really 
the potential for a dynamical scalar field, but rather a nonlinear function of the matter 
variables. This is something of a matter of taste, and we will not pursue this possibility 
in the remainder of the paper. 

We proceed to study structure formation in this model. We start with a Harrison- 



Modified- Source Gravity and Cosmological Structure Formation 



18 




Figure 4. The dimensionless matter power spectrum for MSG and ACDM. As 
expected, in MSG the growth-rate is proportional to k and the power-spectrum 
increases exponentially with k for modes within the horizon. Non-linear scales are 
reached at scales of 240 Mpc (k = 0.004 Mpc -1 ), compared to 10 Mpc= 8/i _1 Mpc 
for ACDM. The linear growth can only be trusted for k < 0.004 Mpc -1 . For smaller 
scales, the theory should return to a GR-like behavior. 

Zel'dovich scale-invariant spectrum of perturbations normalized to COBE by A 2 - = 
5.07 x 10 -5 , (where A is the dimensionless square root of the matter power spectrum 
defined as A 2 = k 3 P(k)/2ir 2 ), and evolve through to today using a numerical code that 
includes both radiation and dark matter. Since we expect small scales to be non-linear, 
it is necessary to find the scale at which linear perturbation theory ceases to be valid. 
These results are shown in figure 0] 

As expected, the non-linear terms are important at much larger scales than in 
ACDM: approximately 240 Mpc compared to 10 Mpc for ACDM. Once the perturbations 
are non-linear, perturbation theory breaks down and an N-body simulation with 
potentials depending on local densities must be used to obtain results. (Nonlinearities 
become important when A ~ 1; the potentials $ and \1/ are safely smaller at that 
time.) However, since for high-enough densities the equations of motion return to GR, 
we expect that the cosmology should behave in a GR-like manner for sufficiently high 
values of 5. 

The evolution of modes at very large scales, for which linear perturbation theory is 
valid throughout the history of the universe, is presented in figures 03 and El Structure 
growth is scale dependent and runaway, as predicted by the approximation in (}55|) . 
with the rate of growth increasing as the universe departs from conventional matter- 
dominated behavior. Such rapid structure formation drives the growth of gravitational 
potentials, which also increase rapidly during the acceleration era in a scale-dependent 
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Figure 5. The evolution of the quantity $ — <J> — observable through the ISW effect — 
for modes close enough to the horizon to remain linear until today. In MSG, this 
combination of potentials grows during dark-energy domination, increasingly rapidly 
for higher k. This would lead to a significant increase in the ISW signal, at least for 
the lowest multipolcs, which are sensitive only to the linear modes at largest scales. In 
LCDM this evolution is scale free. 

manner. This behavior would significantly enhance the Integrated Sachs- Wolfe (ISW) 
effect, at least at the lowest multipoles which are sensitive to only the largest scales, 
which remain linear. In addition, the fact that the gravitational potentials are growing 
at the same time as the density contrasts are increasing would lead to a galaxy-ISW 
correlation of opposite sign to that expected in ACDM; whereas these are correlated for 
ACDM, here they would be anti-correlated. 

7. Conclusions 

We have considered a class of modified-gravity models in which only the constraint 
equation of GR is modified, thereby introducing no new propagating degrees of freedom. 
These "Modified-Source Gravity" models allow for cosmological dynamics in which the 
universe self-accelerates, without the need for dark energy. We demonstrate that there 
exists a class of such theories with a consistent cosmic expansion history and which 
naturally satisfies all solar system tests of gravity. 

As with any explanation for cosmic acceleration, it is important to understand how 
it might be tested and distinguished from other models. To this end, we have studied the 
onset of structure formation in modified-source gravity using linear perturbation theory. 
Fixing parameters so that the background cosmology is well described by MSG, we have 
compared the rate of growth of different Fourier modes of the density perturbations with 
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Figure 6. Exponent of growth rate for comoving density contrast, . In ACDM, it 
is scale- independent once radiation no longer dominates; structure growth slows down 
during dark-energy domination. For MSG, structure growth accelerates. 

those predicted in the ACDM model. We find that, for a given k-mode, growth is more 
rapid in MSG than in ACDM, and therefore that perturbation theory breaks down at 
a correspondingly higher redshift. 

To make detailed progress beyond this point would require an N-body simulation 
with potentials depending on local densities. However, since MSG is constructed so as 
to yield dynamics indistinguishable from GR at high enough densities, we expect this 
rapid growth to cease and to once again resemble that found in ACDM for sufficiently 
high values of 5 m . 

Particularly interesting is the evolution of modes on scales large enough that linear 
perturbation theory is valid throughout the history of the universe. The rapid structure 
formation on these scales drives the growth of gravitational potentials, which increase 
rapidly during the acceleration era in a scale-dependent manner. This behavior enhances 
the Integrated Sachs- Wolfe effect at the lowest multipoles. Since the gravitational 
potentials are growing at the same time as the density contrasts are increasing this 
should lead to a galaxy-ISW anti-correlation, in contrast with that expected in ACDM. 
A natural next step is to attempt to understand the growth of structure in the nonlinear 
regime. 

A primary motivation for this work has been to understand the way in which 
modified gravity can be distinguished from dark energy. Currently, the leading candidate 
for a modified theory of cosmological gravity is the DGP model [21 EH Eh despite lingering 
fundamental issues with the theory Much effort has gone into understanding 

the evolution of cosmological perturbations in DGP gravity, with an emerging consensus 



Modified- Source Gravity and Cosmological Structure Formation 



21 



that the gravitational potentials decay more rapidly at late times in DGP than they 
do in ACDM I2Z1 EHl EH 021 ISHl EH EH EH! In MSG, in contrast, it appears 
as if the potentials generically grow at late times with respect to their conventional 
behavior. It therefore seems to be difficult to imagine a model-independent prediction 
for the way in which modified gravity can be distinguished from theories of dynamical 
dark energy, although simultaneous measurements on the expansion history and the 
evolution of structure do of course provide stringent constraints on any specific model. 
It is clearly important to continue to explore the theoretical consequences of modifying 
general relativity on large scales, to better understand what clues observers should be 
looking for in the quest to solve the puzzle of the accelerating universe. 
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